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Abstract We consider a one-dimensional gas of hard point particles in a finite box that are in ther- 
mal equilibrium and evolving under Hamiltonian dynamics. Tagged particle correlation functions of 
the middle particle are studied. For the special case where all particles have the same mass, we obtain 
analytic results for the velocity auto-correlation function in the short time diffusive regime and the 
long time approach to the saturation value when finite-size effects become relevant. In the case where 
the masses are unequal, numerical simulations indicate sub-diffusive behaviour with mean square dis- 
placement of the tagged particle growing as t/ln(i) with time t. Also various correlation functions, 
involving the velocity and position of the tagged particle, show damped oscillations at long times that 
are absent for the equal mass case. 

Keywords Hamiltonian dynamics • hard particle gas • tagged particle diffusion • velocity autocorre- 
lation function 



1 Introduction 

Observing the dynamics of a single tagged particle in a many particle system constitute a simple way 
of probing the complex dynamics of an interacting many body system and has been studied both 
theoretically 1.2,3,4 and experimentally [SHI]. Much of the theoretical studies on tagged particle 
diffusion have focussed on one-dimensional systems and discussed two situations where the microscopic 
particle dynamics is (i) Hamiltonian QUI] or (ii) stochastic pi [5l[^l[TUl[TTl[T2irT3l . 

For systems with Hamiltonian dynamics, the evolution of the system is completely deterministic 
and all the randomness in the system is due to the randomness in the initial condition. One of the 
earliest result for systems with Hamiltonian evolution is that of Jepsen [1] on tagged particle diffusion 
in a one-dimensional hard particle gas of elastically colliding particles of equal masses. For an infinite 
system at a fixed density of particles Jepsen showed that the mean square deviation (MSD) of a tagged 
particle from its initial position grows linearly with time t. He obtained an explicit expression for the 
diffusion constant and the related velocity autocorrelation function (VAF). This was done by exploiting 
the fact that when two particles of equal mass collide elastically in a one-dimensional system, their 
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velocities are exchanged; if we ignore tags on particles, this is is equivalent to the particles passing 
through each other without colliding, simplifying the dynamics. 

For a finite system, with N particles, there must be corrections to Jepsen's result, since the MSD 
must saturate at long time (to a value that depends on the size of the system). This situation has 
been extensively studied for stochastic dynamics PI WTUIITT1]T2I[T5] but not much for the Hamiltonian 
case (T3UT5TH]. Lebowitz and Sykes [4 considered finite size effects for some special initial conditions. 
In this paper we consider Boltzmann distributed initial coniditions. The first objective of this paper is 
to obtain analytical expressions for the VAF that are valid over the entire regime: both t << N (a la 
Jepsen) and t » N. 

If the particle masses are not all the same in a hard particle gas, there are no analytical results. Since 
the dynamics are expected to be ergodic, the correlation functions should be very different from those 
of the equal mass particle gas. Indeed, a simulation study [16 of a gas where odd and even numbered 
particles have different masses suggested that the decay of the VAF with time t in this model was as 
~ t~ s with 8 < 1 which is completely different from the Jepsen result (~ i -3 ). If this is correct, it 
would imply that tagged particle motion is superdiffusive in this system. The second objective of this 
paper is to accurately obtain the decay of the VAF and other correlation functions for a hard particle 
gas with unequal masses, to see if tagged particle motion is superdiffusive. We perform simulations on 
a one-dimensional gas with alternating masses. To ascertain how robust the numerical results are we 
also do simulations with randomly chosen masses. 

Although there has been considerable work on the (hydro)dynamics of one dimensional hard particle 
gas and other systems in the context of heat conduction [17] , this involves the propagation of conserved 
quantities as a function of position and time without reference to the identity of each particle. This 
changes things considerably: for instance, conserved quantities propagate ballistically for an equal 
mass hard particle gas, resulting in a thermal conductivity proportional to N, while tagged particle 
dynamics in the same system is diffusive. Thus here we approach the dynamics from a perspective that 
is different from the heat conduction literature. 

In sec. ^ we define the model and dynamics and give analytic results for the VAF in the special 
case where all masses are equal. In sec. ^ we present the simulation results for the VAF and other 
correlation functions for the general case where masses are not all equal. We summarize our results in 
Sec. |4| Some details of the calculation are given in Appendix [Aj 



2 Analytic results for equal mass hard-particle gas 

Here we consider a gas of N = 2M + 1 point particles in a one-dimensional box of length L. The 
particles interact with each other through hard collisions conserving energy and momentum. The 
Hamiltonian of the system thus consists of only kinetic energy. All the particles have the same mass m. 
In any interparticle collision, the two colliding particles exchange velocities. When a terminal particle 
collides with the adjacent wall, its velocity is reversed. The initial state of the system is drawn from 
the canonical ensemble at temperature T. Therefore, the initial positions of the particles are uniformly 
distributed in the box. Let Xi be the position of the i-th particle measured with respect to the "left" 
wall, and < x\ < X2 < ■ ■ ■ < %n-i < x n < L. The initial velocities of the particles are choosen 
independently from the Gaussian distribution with zero mean and a variance v 2 = ksT /m. 

By exchanging the identities of the particles emerging from collisions, one can effectively treat the 
system as non-interacting pQ. In the non- interacting picture, each particle executes an independent 
motion. The particles pass through each other when they 'collide' and reflect off the walls at x = 
and x — L. The initial condition is that each particle is independently chosen from the single particle 
distribution p(x,v) = L- 1 {2i:v 2 )- 1 / 2 e~ v ' ' l 2 ^ , where v 2 = k B T/m. To find the VAF of the middle 
particle in the interacting-system from the dynamics of the non-interacting system, we note that there 
are two possibilities in the non-interacting picture: (1) the same particle is the middle particle at 
both times t = and t, or (2) two different particles are at the middle position at times t = and t 
respectively. We denote the VAF corresponding to these two cases by (v m (0)v m (t)) l and (um(0)"mW}2 
respectively. The complete VAF is given by (vm(O)vm(^)) = (vm(0)vm(£))i + (um(0)umW)2- We now 
present a physically motivated derivation of these two quantities. A direct derivation and some more 
details are given in an appendix. 
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We first define a few quantities. The probability density for a (non-interacting) particle to be at x 
at time t = and y at time t is 



P(x,y;t) = 
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The first line is easily obtained by realizing that, for a free particle in an infinite box with a Gaussian 
velocity distribution, the corresponding probability density is (Lvty/Zn^ 1 exp[— (x — y) 2 /(2v 2 t 2 )], and 
the boundaries at x = and L set up an infinite sequence of image sources. The second line is obtained 
using the Poisson resummation formula or by realizing that with r = i 2 , the first expression satisfies 
d T P = —v 2 d 2 P with initial condition P(x, y, 0) = S(x — y)/L and boundary conditions d y P(x,0,T) = 
d y P(x, L, t) = 0, and expressing this in terms of the eigenfunctions of the Laplacian. As a variant of 
Eq. (fil) we also define the function 
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In terms of these functions, the correlation function («m(0)vm(£)}i can be found by picking one 
of the non-interacting particles at random, calculating the probability that it goes from (x,Q) to (y, t) 
and that it is in the middle at both t — and t, multiplying by v(0)v(t) and integrating over x and y. 
The multiplication by v(0)v(t) is equivalent to inserting a factor of (2nL + y — x) 2 /t 2 in the first term of 
the first line of Eq.([lJ) and a factor of —(2nL — y — x) 2 /t 2 in the second term, since they correspond to 
even and odd number of reflections respectively. Thus one obtains the normalized correlation function 
(see appendix) 



VV 
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where P^ (x, y, t) is the probability that there are an equal number of particles to the left and right 
of x and y at t — and t respectively. 

Turning to («jvf(0)«Af(i))2 ) we pick two particles at random at time t = 0, calculate the probability 
that they go from (x, 0) to (y, t) and (x, 0) to (y, t), that there are an equal number of particles on both 
sides of x and y at t = and t respectively, multiplying by v(0)v(t) and integrating with respect to 
x, y, x, y. From Eq.([T]), multiplying P(x, y, t) by u(0) is equivalent to inserting a factor of (2nL + y— x) jt 
and (2nL — y — x)/t in front of the first and second terms respectively in the first line. Also, multiplying 
P(x, y, t) by v(t) is equivalent to inserting a factor of (2nL+y — x)/t and —{2nL — y — x)/t in front of the 
two terms. Converting these factors to appropriate derivatives, we have for the normalized correlation 
function 



C® = N{N-1) 



dxdxdydy P^ (x,x,y,y, t) [vt d x P(x, y, t)] [-vt dyP(x, y)] . 



(4) 



(2) 

where P^ (x,y,x,y,t) is the probability that there are an equal number of particles on both sides of 
x and y at t = and t respectively, given that there is one particle at (x, 0) and at (y,t). 

(1 2) 

To proceed further, we need the expressions for P^' . For this we define p \.(x,y;t) as the prob- 
ability that a particle is to the left of x at t — and to the right of y at time t. Let , p and p++ 



4 



be similarly defined. Thus 



dx' / dy'P(x',y';t) , 
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In terms of the expressions defined in Eqs.([5|, it is straightforward to see that 
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where we used the identity p++ +p~\ +p \- + p = 1, the fact that N — 1 is even and the integrands 

are unchanged if 9,4> are increased by tt. The angular integrals enforce the conditions that if m + n 
particles are to right of x at time t = 0, of which n particles cross from right to left in time t, then n 
particles cross from left to right and m particles remain on the left, so that the number of particles on 
both sides of x at time t — and y at time t is m + n (see appendix for more details). From Eqs.((3| 
and ([6| we get 
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Using similar arguments as used for Eq.(|6| (see appendix), one can write 
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Using the second line of Eq.([T]) for P(x,y,t) and P(x,y,t), integrating over y and x, and comparing 
to the second line of Eq. ^ we obtain 
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2.1 Short time regime 

When vt << L, the tagged particle does not feel the effect of the walls and we can make the following 
approximations 



P(x,y,t) = P-(x,y,t) 
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In this limit the expressions for p [_, etc. given in Eq. Q also simplify by using Eq. (11) and taking 

the limits of the y' integral to be from y to oo for p |_,p++ and from — oo to y for p_| , p We then 

get 
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Hence we get iV(p+ + +p__) = N-2Xq{z_), = -Xz+, N{p + -+p_ + ) = 2Ag(z_), N(p + _- 

p |_) = — Az_ . Using these in Eqs. Q and ( 10 ), and changing variables from xo, Xt and xq, it to z + , z_ 

we get, for large N 



cilHt) 



cilHt) 



A 



/ 2^{2tt) 2 

x e" 
A 2 '°° 



dz4 



/■tt/2 ptt 

dz-z 2 _e- z - /2 / d<t> / d6 

-tt/2 J— it 



oo 

N( l — cos (fx) —iXzjf. sin <f> —2\q{z — ) (1 — cos 9) — i\z— sin £ 



4.3 ' dz+ 



— oo J — oo 



dz-e 



tt/2 



tt/2 J-tt 



d9 (cos 4> — cos 0) 



v — iV(l— cos 0) — i\z+ sin </> p — 2Ag(js_ ) (1 — cos $) — iXz— sin t 



For large N, the major contribution of the integral over <f> comes from the region around <f> = 0. 
Therefore, the <f> integral can be performed by expanding around <j> = to make it a Gaussian integral 
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Fig. 1 Plot of the separate contributions and CvV to the velocity-autocorrelation function of equal mass 
hard-particle gas, for two different system sizes (N = 101, 201) with fixed density p = 1, and v = 1. The solid 
lines cor respond to the simulation data, whereas the points are from the analytical results given by Eqs. |7| 
and (10 1 for short time behaviors (a), and Eq. (261 for the long time behaviors (■). In the first panel, tne 
dashed rectangle is enlarged in the inset. 



(while extending the limits to ±oo). Subsequently, one can also perform the Gaussian integral over z + . 
This leads to the following expressions: 
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Thus we have closed form expressions of the VAF which are valid in the entire short time regime. For 
any value of A = pvt, these integrals can be performed numerically, and as we see from Fig. [T] the 
results are in excellent agreement with the numerical simulation. 

We now analyze the above expression in the large A limit, i.e., (vp)" 1 -C t. Together with the 
condition vt << L for being in the short time regime, this means that t is much larger than the typical 
time between interparticle collisions, i.e. outside the ballistic regime, while being much smaller than 
the time it takes to see finite size effects. We first make a change of variables \f\z^ = z and \f\6 = x. 
The integrands can then be expanded as a power series in powers of 1/A. The integrals acquire the 
forms: 

C vv{ t ) = 77TW2 dz dxe'^" lxz J2a n (x,z)X- n 1 (16) 

[Z7T) J-oo J - VXtt n=2 
-i poo py/X'K 2 °° 
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where a n (x,z) and b n (x, z) are polynomials in x and z. For example, a,2(x, z) = z 2 , 61 (x, z) — x 2 /2, 
and so on. Now, integrating term by term (while extending the integrating limits of x to ±00) we get 

^ ) W = ^A- 2 -fc|A- 3 + 0(A- 5 ) ) (18) 

C$(*) = -^A- 2 - A- 3 + 0(A- 4 ). (19) 

Therefore, adding the above two results, we recover Jepsen's result pQ 

a„(i) = -|^A- 3 + 0(A- 4 ). (20) 



2.2 Long time regime 

In the limit vt » L we integrate Eqs.([5]) using the second line of Eq.Q: 
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Expanding around x = y = L/2 we get to leading order N(p ++ + p ) = N[l/2 + 2a(t)], N(p ++ — 
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Using these and the expression of P_ from Eq. ^ in Eqs. Q and (10), we find the following results 
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Performing the Gaussian integrals we find that C\, v vanishes and hence to 0(1/N) the velocity auto- 
correlation is given by 

fc=l,3,5.. ' 



As seen from Fig. [T] the above expression describes the numerical simulation data very well. The late 
time behaviour C vv (t) ~ exp(— ir 2 v 2 t 2 /2L 2 ) was earlier obtained in [13] . 



3 Simulation results 

As mentioned earlier, there are no analytical results when the particle masses in the one dimensional 
gas are not all equal. We turn to numerical simulations for such systems; the simulations also confirm 
the analytical results of the previous section, as shown in Figure [TJ The Hamilt onian for the system is 
H = y^i—i ^mi±f, with < X\ < x 2 ■ ■ ■ < < L. After an elastic collision between two neighboring 
particles (say I and I + 1) with velocities Vi, vi + i and masses mi, mi + i respectively, they emerge with 
new velocities v\ and v' l+1 . From momentum and energy conservation we have: 
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Between collisions the particles move with constant velocity. 

We simulate this system using an event-driven algorithm and compute the correlation functions 
([Ax(t)] 2 }, (Ax(t)v(0)) and (v(t)v(0) of the central particle, where Ax(t) = xm{t) — xm (0) and v(t) = 
VM(t)- The average (• • •) is taken over initial configurations chosen from the equilibrium distribution, 
where the particles are uniformly distributed in the box with density p — N/L, while the velocity of 
each particle is independently chosen from the distribution (m/27rfcsT) 1 / 2 e~ mv / 2fcfiT . Note that the 
three correlation functions are related to each other as 



t 

{[Ax{t)f) = (Ax(t)v(t)) = (Ax(t)v(0)} = I (v(Q)v(t'))dt' = D{t) . 



2dt 



When the tagged particle shows diffusive behaviour then lim t _ J . 00 D(t) reaches a constant value for an 
infinite system and this gives the diffusion constant. On the other hand for sub-diffusion D(t) vanishes 
as t — > oo whereas for super-diffusion it diverges. 

Just as for the equal mass system, for any finite system of size L there is a short time regime 
during which the tagged particle at the centre does not feel the effect of the boundaries and during 
this time, correlation functions have the same behaviour as the infinite system. The time at which 
the system size effects start showing up is given by t sat ~ L/c s , where c s — y/3P/ ' p m is the adiabatic 
sound velocity in the hard particle gas, with P the pressure and p m the average mass density. For our 
numerical simulations, P = pksT = 1 and p m = 1, which gives c s = a/3- We now present the results 
for the correlation functions in the short-time and long-time regimes. 

In Fig. [2j we show the simulation results for the correlation functions for a one dimensional hard 
particle gas with masses that alternate between 1.5 and 0.5. Here the data is shown for the case where 
the tagged particle has mass 1.5, and similar results are obtained for the case when the tagged particle 
is lighter. For comparison, the results for an equal mass gas are also shown. After the expected initial 
ballistic regime, the MSD (Ax 2 (t)) grows approximately linearly, indicating roughly diffusive motion. 
Simulation results of Marro and Masoliver [TB] obtained (v(0)v(t)) ~ — t~ s with 8 < 1 for the gas 
with alternating masses, which would imply (slightly) superdiffusive behavior. It is easiest to notice 
any deviations from diffusive behavior in the plot of (Ax(t)v(0)) , where diffusive or superdiffusive 
behavior would correspond to (after the ballistic regime) a horizontal or rising straight line respec- 
tively. Instead, Figure [2] shows that (Ax(t)v(0)) decreases as t is increased beyond the ballistic regime, 
implying subdiffusive behavior. This is seen more clearly in Fig. ^ where we observe the dependence 
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Fig. 2 (color online) Various correlation functions for alternate mass hard particle gas (solid lines) with 
N = 101 (blue) and N = 201 (red) particles, density p = 1 and ksT = 1. The alternate particles have masses 
1.5 and 0.5 and in this simulation, the middle particle had mass 1.5. The data is obtained by averaging over 
10 9 equilibrium initial conditions. For comparison, the correlation functions for an equal mass gas with masses 
1 is also shown (dotted lines). 



(Ax(t)v(O)) ~ a/(b + hit). This would correspond to an MSD whose leading part in this regime is 

- t/hxt and a VAF l/(tln 2 t). As seen in Fig. ^ the logarithmic corrections are difficult to 

observe in the MSD and VAF. The smallness of the deviation from diffusive behavior implies that the 
apparently linear dependence of (Ax(t)v(0)) on lni, observed in [16] ■ is nevertheless consistent with 
the above observed form for small Int. However we note that the linear logarithmic dependence on 
time as proposed in [16] cannot be valid at large times - and therefore our form is more appropriate. 
Other functional forms for are also possible. For instance, (Ax(t)v(O)) ~ t~ a with a very small a would 
also imply the subdiffusive behaviour (Ax 2 (t)) ~ t x ~ a . 

At long times, the effect of finite size of the box sets in and the MSD saturates: ([x(t) — a;(0)] 2 ) — > 
2([x(t) — L/2] 2 ), which can be easily evaluated in equilibrium to be L 2 /(2N), independent of the 
particle masses in the gas. We observe this in Fig. The main difference between the equal mass and 
alternate mass systems is that the MSD for the equal mass case approaches its saturation value without 
oscillations, while for alternate mass case there are damped oscillations as saturation is approached, 
while always remaining below the MSD for the equal mass case. The oscillations in the alternate mass 
system also show up in the other two correlation functions. 

The oscillations in the MSD are seen more clearly in Fig. |1J where the data is plotted differently. 
The period of oscillation is proportional to N, in agreement with our discussion earlier in this section 
where they were ascribed to sound waves reflecting from the boundary, which takes a time ~ L/c s . 
However, the amplitude of the oscillations does not show a simple scaling with N; it is clear from the 
figure that they are damped out in fewer cycles for smaller N, making it impossible to collapse the 
data onto a single curve by rescaling the vertical axis. 

To check for the robustness of our results we have also performed simulations of a gas with random 
distribution of masses. Each particle was assigned a mass from a uniform distribution between 0.5 and 
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t 



Fig. 3 Plot of D(t) = (Ax(t)v(0)) for the alternate mass gas for various system sizes. We see clearly the 
logarithmic decay of the diffusion constant. The parameters for the shown funciton are a — 8 and b = 19. 
For comparision we also show the corresponding equal mass data (dashed line) which shows saturation to the 

expected Jepsen value l/v27r ~ 0.4. 




1.5. We looked at tagged-particle correlations of the central particle whose mass was fixed at 0.5. The 
correlations fluctuate between different mass realizations and we took an average over 32 realizations. 
The results are plotted in Fig. ([5| where we see the same qualitative features as for the alternate mass 
case. 
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Fig. 5 Various correlation functions (in the short-time regime) for random mass hard particle gas with 
N = 101 and N = 201 particles and density p = 1. The mass of the middle particle was always taken to be 0.5 
and the results are an average over 32 different random mass realizations. 



4 Summary 



We have studied tagged particle correlations of the middle particle in a system of N hard point particles 
confined in a one-dimensional box of length L and in thermal equilibrium. For the case where the 
masses of all particles are equal we obtained analytic results for the finite-size velocity auto-correlation 
function using the approach of Jepsen. We have presented a somewhat simpler and physically motivated 
calculation of the velocity auto-correlation function and obtained closed form expressions valid at both 
short times (including the ballistic and diffusive regimes) and long times (when finite size effects show 
up). While here we have only presented results for the velocity auto-correlations, it is straightforward 
to obtain other correlation functions using our approach. 

Next we have presented simulation results for the case of a hard-point gas where the particles 
have unequal masses. Two cases are studied, one where particles have alternate masses and the other 
where the masses are random. In both cases we find that the behaviour of correlation functions is 
qualitatively different from the equal mass case. The correlation (Ax(t)v(0)) does not saturate to a 
constant (expected for a diffusive behaviour) and instead shows a slow decay consistent with the form 
(Ax 2 (t)) ~ t/hxt. Correspondingly the VAF decays as ~ l/(iln 2 i) which is completely different from 
the equal mass form ~ l/t 3 . This indicates that tagged-particle motion is sub-diffusive. However it 
is difficult to see this sub-diffusive behaviour directly in the mean square displacement of the tagged 
particle since the deviation from linear time-dependence is small. These results are surprising since 
simulations with other interacting systems such as Lennard- Jones gases have found diffusive motion 
and l/t 3 decay of the velocity auto-correlation function [18]. Understanding this difference as well as 
studying tagged particle motion in other interacting systems and higher dimensional systems remain 
interesting open problems. 
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A Details of calculation 

In this appendix, we provide a more detailed calculation of the velocity autocorrelation function for the hard 
particle gas of equal mass. This is an alternative to the derivation of some of the key equations in this paper. 
To compute (om(0)dmM)i, we pick at time t = one of the non-interacting particles at random from the 
distribution p(xo , Do) ■ At time t let the position and velocity of the particle be given by xt (xo , vo ) and Vt (xo , vo ) 

respectively. We then calculate the probability, P^\xo, x t ), that it has an equal number of particles to its left 
and right at both the initial and final times, i.e., at t — and t. For {vM(0)vM(t))2, we pick two non-interacting 
particles from the distribution p(xo,Vo) p(xo,Vo) and let them evolve to x t ,v t and Xt,v t respectively. We then 

calculate the probability Pju (xo, it) that at time t = 0, the first particle xo is the middle particle while at 
time t the second particle x(t) is the middle particle. The normalized VAF is thus given by: 



V 2 V 2 Jo L J.^ V27TU 



C v J (t) 



(VM(t)VM(0)}2 



f dxo 
J ~L 



dxo dvodvo 
T~ 2-irv 2 



V V t e -("0 2 +^)/2- 2 pW ( XOtXt ,x ,X t ,t) 



(28) 



(29) 



These forms together with the explicit expressions of P^'^ di scus sed below, agree with those given in [4j. 



We now make a change of variables from xo,vo to xo,Xt in Eq. (281 and from xo,vo,Sto,vo to xo, Xt, xo, Xt in 
Eq. (29}. 

In the non-interacting picture, xt and vt, as well as the number of collisions m, suffered by the particle with 
the walls upto time t, are completely determined by the initial configuration (xo,Vo)- The number of collisions 
with the wall is given by 

' xo + v t 

if v > 0, 



L 



L — xo — vot 
L 



(30) 



if v < , 



where [^J is the integral part of x. When m is even, we have vt = vo whereas Vt = —vq for odd m. The 
final position x t is given by one of the following relations depending on m and vq. When m is even, we have 
xo + vot = mL + xt for no > and L — xo — vot = mL + L — xt for no < 0. On the other hand for odd m we get 
xo + vot = mL + L — xt for vo > and L — xo — vot — mL + x t for vo < 0. Combining all these four cases, we 
can write xo + vot — 2nL ± xt- Here n = m/2 and — m/2 respectively for the first two cases where m is even 
and the plus sign is taken. For the last two cases, where m is odd, n = (m + l)/2 and — (m + l)/2 respectively 
and the minus sign is taken. In other words, for a given values of xo and vo in the relations xo + vot — 2nL ± x t 
and vt = ±«o, the values of n and xt, and the signs taken from the ± are uniquely determined. Therefore, 
inse rting the term [S(xo + vot — 2nL — xt)5(vt — vo) + S(xo + vot — 2nL + xt)8(v t + «o)] in the integrand of 
Eq. ( 28 1 while integrating over x t and vt, and summing over all integer values of n, does not change the result, 
i.e., 



(v M (t)v u mi=N f dx t r dv t V f ^ I" ^e-"2 /s,J W F 1 ? ) (io,x t ,t) 
Jo J-oo n tz aa Jo L J-oo V2nv 

x [S(x + v t — 2nL — x t )3(v t — v ) + S(x + v t — 2nL + x t )S(v t + v )] . (31) 
Now, carrying out the integrations over v t and vo, after some straightforward manipulation we obtain 

C^^n( dx [ dx t P£\xo,xt,t)v d„ P-(x ,x t ,t), . (32) 
Jo Jo 



The second part of the velocity autocorrelation function is given by Eq. (29 1 and in this c ase we trade the 
vo,i>o integrals for xt, x t by introducing two sets of 5- function, one for each particle as in Eq ( |31[ ). After some 
manipulations we then get 

C$=N(N — 1) J "J dxodx dx t dxt P^\xo,x ,xt,x t ) 

[vt d X0 P(x ,x t )] [-vt di t P(x ,x t )], (33) 

Evaluation of P^ (xo, x t ): This gives the probability that, at t = and at time t, the selected particle has 
an equal number of particles to its left and right. We note that the remaining N — 1 particles are independent 
of each other and the selected particle. Let p — ^(xo,x t ;t) be the probability that one of these particles is to 
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the left of xo at t = and to the right of x t at time t. Let p^ , p and p++ be similarly defined. In terms of 

these probabilities, it is easily seen that 



2^ n 1 !n 2 !n 3 !n 4 ! ;P — P - +P+ - P++ ni '" 4 ™ 2 ' n3 ' 1 j 

n 1 +7i2+n 3 +n 4 — N — l 



where in the summand, m particles go from the left of xq to the left of x t , ni particles from the left to the 
right, ri3 particles from the right to the left, and 774 particles from the right to the right. The two Kronecker 
delta functions ensure that an equal number of particles cross the selected particle in both directions in 
time t and that an equal number of particles remain on either side of the selected particle. Together, these 
conditions are equivalent to an equal number of particles being on either side of the selected particle at time 
and t, that is, 77 1 + n% = 773 + 714 and 711+713 = 772 + 774. The multinomial coefficient takes care of all 
possible permutations among the particles. Now, using the integral representation of the Kronecker delta, 
S m ,n = (2tt) _1 J Q 2Tr e l ( m - n )° fiQ m the above equation immediately gives Eq. Q 

Evaluation of P^\xo, x t ): In calculating pffi we have to keep track of both the particles. There arise 
four situations: (a) xo > xo and xt < xt, (b) xo > xo and xt > xt, (c) xo < xo and Xt > xt, and (d) xo < xo 
and x t < xt- Let there be m particles go from the left of Xo to the left of xt, ni particles from the left 
to the right, 773 particles from the right to the left, and 774 particles from the right to the right. Since two 
of the particles are considered separately, the rest can be chosen (N — 2)!/(7ii!n2!773!n4!) different ways and 
7ii+772 + 773 +7i4 = N— 2. Now, in the first situation we have (a) m +712 + 1 = 773 + 774 and 771 +773 + 1 = 772+774. 
These conditions are equivalent to 712 = 774 and 77 1 =774 — 1. Similarly one can work out the conditions for the 
other three situations which gives (b) 771 = 774 and 712 = 773 — 1, (c) 772 = 773 and 771 = 774 + 1, and (d) n\ = 774 

and 77 2 = 773 + 1, respectively. Following the procedure used to evaluate P^\ we can easily find pjp as given 
by Eq. pj, where the extra phase factor tp(0, <j>) originates from addend ±1 that appear in the relations among 
ni's above, and ijj(d,<j)) = e _l<#> , e~ %e , e** and e l6 for situations (a), (b), (c) and (d) respectively. 

Evaluation of P(xo,x t ): The joint probability density function for a (non-interacting) particle to be 
between x and x + dx at t = and between y and y + dy at time t is given by 

P{x,y) = {5(x - x )5(y - x t )) 

rL dx f 00 e - v o/ 2 v 2 

—— / dv — , — _ S(x — x )S(y — x t ) 
Li i_oo \/2nv 

= f dx t [ dv t y2 [ f dv 6 7 — _ S(x - x )S(y - x t ) 

JO J -co n T^ oc J0 L J -co V27T7; 

x [5(xq + v t — 2nL — x t ) S(v t — v ) + 5(x + v t — 2nL + x t ) S(v t + v )] . (35) 
Now, carrying out the integrations over all the variables gives the first line of Eq. |l]). 
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